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Abstract 

This paper presents buckling and postbuckling re- 
sults for simply supported aluminum plates loaded 
in compression. Buckling results have been plot- 
ted to show the effects of thickness on the buckling 
stress coefficient. Buckling results are given for var- 
ious length-to-width ratios. Postbuckling results for 
plates with transverse shearing flexibility arc com- 
pared with results from classical theory for various 
width-to-thickness ratios. The plates are considered 
to be long with side edges simply supported and 
free of stress. The plates are subjected to a lon- 
gitudinal compressive end-shortening displacement. 
Characteristic curves indicating the average longi- 
tudinal direct stress resultant as a function of the 
applied displacements are calculated based on four 
different theories: classical Von Kannan theory, first- 
order shear deformation theory, higher order shear 
deformation theory, and three-dimensional flexibil- 
ity theory. Present results indicate that the three- 
dimensional flexibility theory gives the lowest buck- 
ling loads and, therefore, the most accurate results. 
The higher order shear deformation theory has fewer 
unknowns than the three-dimensional flexibility the- 
ory but is not as accurate. The* figures presented 
show that for postbuckling small differences occur 
in the average longitudinal direct stress resultant, in 
the maximum values of the transverse stress resul- 
tant and shear stress resultant, and in the maximum 
transverse displacements calculated when the effects 
of transverse shear flexibility are included in the var- 
ious theories. 

Introduction 

The increasing interest in minimum weight de- 
signs for aeronautical and aerospace structures has 
generated substantial interest in the analysis of the 
elastic stability and postbuckling behavior of struc- 
tures subjected to compressive in-plane loads. For 
thin homogeneous plates, classical plate theory pre- 
dicts deformations and in- plane stresses that are 
comparable to those of three-dimensional elasticity. 
Transverse stresses in thin plates are generally small 
compared to in-plane stresses, and thus both clas- 
sical theory and first-order shear deformation the- 
ory give satisfactory results. However, since both 
theories are two-dimensional, they are not accurate 
enough to predict transverse stresses directly. Accu- 
rate nonlinear theories are required for the analysis 
of thick plates, in which these transverse stresses be- 
come more significant. 

It is often sufficient to use an accurate non- 
linear two-dimensional theory to solve some three- 
dimensional nonlinear elasticity problems. One such 


theory has been derived in reference 1 for laminated 
and thick plates with three-dimensional flexibility ef- 
fects. This theory can predict directly the trans- 
verse st resses as well as the in-plane stresses by using 
trigonometric terms in addition to the usual constant 
and linear terms representing through-the-thickness 
variation of the displacements. However, this theory 
cannot satisfy the surface boundary conditions of a 
plate. 

The purpose of the present paper is to present the 
results of an investigation of the buckling and post- 
buckling response of isotropic plates loaded in com- 
pression using classical nonlinear Von Karman theory 
and thiee nonlinear transverse shearing theories and 
to compare results for these four theories for differ- 
ent values of plate width-to-thickness ratios in the 
postbuckling range. The nonlinear transverse shear- 
ing theories are first-order shear deformation theory 
(refs. 2 and 3), higher order shear deformation the- 
ory (ref. 4), and three-dimensional flexibility theory 
(ref. T . The present derivation of the higher or- 
der shear deformation theory has the advantage of 
having nonlinear through-the-thickness terms with- 
out contributing additional unknowns to first-order 
shear deformation theory. In addition, it satisfies the 
surface boundary conditions of the plate. The es- 
sential difference between the theories is the use of 
cubic or trigonometric terms in addition to the con- 
stant and linear terms that represent the through- 
the-thickness variation of the in-plane displacements. 
The present paper presents the derivation of the non- 
linear plate equations for buckling of plates loaded 
in axial compression for both higher order theories. 
This paper also presents postbuckling results for the 
average longitudinal compressive direct stress resul- 
tant and maximum stress resultants jus a function 
of the applied displacements, and maximum out-of- 
plane displacement as a function of the applied dis- 
placement. The plates considered in this paper are 
long with side edges simply supported, free of nor- 
mal stress, and free to slide along the edges to give 
constant strain. Results of the four theories arc for 
aluminum plates loaded in uniaxial compressive end 
shortening. 
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plate extension al 
stiffnesses 

dimension of rectangular 
plate parallel to X,Y , 
and Z axis, respectively 

stiffnesses used in 
Hooke’s law 
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plate bending stiffnesses 

E 

Young’s modulus 
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plate stiffness components 
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Lxi /> Myy M X y 

moment resultants in the 
plate 
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in-plane stress resultants 
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Nyz 1 N X z 

transverse stress resul- 
tants in the plate 
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average axial stress 
resultant 

N X .CT 

value of N X SiV at buckling 
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functions of y defined in 
appendix 
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applied end-shortening 

u,v, w 

displacement in x , and 

z direction, respectively 

x . y , 2 

plate coordinates 

0\ 1 ft 2 

functions of x and y 
defined in the appendix 

f x 5 ( y * f z ’ 

7 yz^ Ixzi ~fxy 

strains in the plate 

A 

buckle half- wavelength 

V 

Poisson’s ratio 

Vyi &Z1 
T yz 1 Txz f Txy 

stresses in the plate 

Theory 



A brief outline of the derivation of the four differ- 
ent theories compared in this paper is presented in 
this section. The derivation of equations using clas- 
sical Kirchhoff theory has been presented in refer- 
ence 5 . The derivation of equations using first-order 
shear deformation theory has been presented in refer- 
ences 2 and 3 . The two higher order theories are not 
given in detail elsewhere, so they are presented in the 
appendix. The general. approach used in deriving the 
equations to be solved is the same as in reference 6. 
First, the displacement functions for each theory are 
identified. The nonlinear strain-displacement rela- 
tions and the assumption that the displacements are 
sinusoidally periodic along the length of the infinitely 
long plate are incorporated. Stress-strain relations 
are defined for a “specially orthotropic” plate, and 
application of the principle of virtual work leads to 


ordinary differential equations and variational lv con- 
sistent boundary conditions, which are solved by a 
computer program based on Newton’s method. 

The displacements considered for each theory arc 

Classical Kirchhoff theory: 


u(x, y, z) = u°(x, y) - w° (2, y) \ ' 

h 

v{x,y,z)=v°(x,y)-w°(x,y)^ > (!) 

w(x, y, z) = w°(x, y) 

First-order shear deformation theory: 


u(x, y , z) = u°(x, y) + u a (x , y) \ ' 

h 

v{x,y,z) = v°(x,y) +v a {x,y)y > 

h 

w(x, y, z) = w°(x, y) 

Higher order shear deformation theory: 


( 2 ) 


u{x,y,z) = u°(x, y) + { u a (x,y) - - 


|«°( 


Ah 


u a (x, y) 


+ w.° (x.y) 


ay 


v{x,y,z) = 1 J°(x,y) + { v n (x,y) 


z 

' h 

Ah 


v a (x.y) 


( 3 ) 


+ W -y (Z-y) 
w{x,y,z) = W u (x,y) 


(I)’ s 


Three-dimensional flexibility theory 

u(x, y , z) = u°(x, y) A- u"(x, y) ~ + u s (x, y) sin ~ 'j 


v(x, V , z) = v°(x, y) + v“(x, y) ^ + V s (x. y) sin 7 

n h 


) ( 4 ) 


w(x,y,z) = w°(x,y) + w r (x,y)< 


Both the classical Kirchhoff and the first-order shear 
deformation theories have in-plane deformations u 
and v which are linear in 2. Classical theory, however, 
has the additional assumption that there is zero 
transverse shearing (7 xz = ^ yz — 0), thus eliminating 
u a and v a in favor of derivatives of w ° . 
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Higher order shear deformation theory considers 
in-plane deformations u and v which are cubic in 0 . 
As explained in reference 6, the squared- in - 2 term 
vanishes and the cubic term does not introduce any 
new variables beyond those that appear in first-order 
shear deformation theory if the boundary conditions 
are satisfied at 2 = ±hf 2 . The three-dimensional 
flexibility theory considers trigonometric terms in 
u, u, and w beyond the expressions considered for the 
deformations of first-order shear deformation theory. 
In this paper, the superscript o corresponds to the 
constant-in-2 terms, the superscript a corresponds to 
the algebraic-in -2 terms, and the superscripts s and 
c correspond to the trigonometric-in-2 terms. 

To account for the applied displacement U 


u°{x,y ) = -U^ +Ul{x,y) 
v°(x,y) = v%(y) + V'?(x,y) 


(5) 


where numbered subscripts for loads, displacements, 
and curvatures indicate a y-depcndence only. To 
satisfy the assumption that the displacements are 
sinusoidally periodic along the length x 


Hooke's law relating stresses to strains for a “spe- 
cially orthotropic” material is used here 


°x 


~C n 

C12 

0 

0 

0 

0 


( X 

a V 


Cn 

Q22 

0 

0 

0 

0 


( y 

<?Z 


0 

0 

C33 

0 

0 

0 


t z 

Tyz 


0 

0 

0 

C44 

i) 

0 



~xz 


0 

0 

0 

0 

C55 

0 


Ixz 



0 

0 

0 

0 

0 





( 8 ) 


Ordinary differential equations and variationallv 
consistent boundary conditions are derived using the 
principle of virtual work, and the equations are solved 
by Newton’s method in a computer program dis- 
cussed in reference 7. The principle of virtual work 
applied to the internal forces of a three-dimensional 
body considered here is 

ra rb fh/2 

6U= / / (ff xfax + fey + cr 2 fez 

Jo Jo J-h/2 

T Ty Z b-fyz + T XZ b^ XZ T Try b^y) Az dlj Ax 

0) 

and t .he set of simple support boundary conditions at 
y — 0 and y — b used are 


2 :T.r ^ 


U2 — r*2 iv) s* 11 - 


n n/ x 2 KX . 

v 2 =v$(y)cos--~ J 


( 6 ) 


u = U ' 

*y = 0 > 
w = 0 

My= 0 . 


( 10 ) 


All the other u coefficients can be expressed as func- 
tions of y multiplied by cos 7r^/A, where A is the 
half-wavelength of the buckled plate. All the other 
v and w coefficients can be expressed as functions of 
y multiplied by sin nx/X. The strain-displacement 
relations used are 


The half- wavelength A of the assumed deformations 
for the long plates considered is chosen to minimize 
the buckling load for each given applied deformation. 

Results and Discussion 


1 2 

+ 2 W ' X 

1 2 

e y = v iy + 2 w ’y 
€ Z = W, z 
lyz — V)Z 
Ixz = TtC,x 

= +v,x 


> 


w 


y ) 


The results obtained in this study for aluminum 
plates with geometry as shown in figure 1 are based 
on the mechanical properties of Young's modulus 
E = 10.7 x 10 6 psi and Poisson’s ratio y = 0.33. 
Buckling results presented in figure 2 for an alu- 
minum plate show the variation of buckling stress 
coefficient with width-to-thickness ratio b/h given by 
the four theories for a range of length-to-width ra- 
tio a/b. The differences in the buckling results for 
aluminum plates with width-to-thickness ratios less 
than 10 illustrate the need to include the effects of 
transverse shearing in determining the compressive 
buckling stress. 
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Characteristic curves for different values of b/h 
are presented in figure 3 for the postbuckling re- 
sponse of rectangular aluminum plates loaded in 
compression. For a b/h value of 100, only one curve 
is shown, since the corresponding results for each 
theory are approximately the same. For lower val- 
ues of b/h, figure 3 shows that the normalized end- 
shortening U is larger for a shear deformable plate 
than for its infinitely rigid counterpart for a fixed 
value of normalized stress resultant N x/av . Present 
results indicate that three-dimensional flexibility the- 
ory gives lower buckling loads and therefore more ac- 
curate results. The shear deformation theory results 
arc converged solutions for all values of b/h inves- 
tigated, and this theory has the advantage of fewer 
unknowns than the three-dimensional theory. Re- 
sults for the normalized stress resultants N yj rlax and 
iVj^max) for different values of b/h, are presented 
in figures 4 and 5, respectively, and corresponding 
results for the normalized maximum deflection w max 
are presented in figure 6. The results presented in fig- 
ure 4 show that the higher order shear deformation 
theory gives the lowest value of normalized N y iniiX 
for a given value of normalized end-shortening U . 
The results presented in figure 5 show that the 
normalized N xy>uvdx as a function of the normalized 
end-shortening U is nearly independent of the width- 
to-thickness ratio b/h for the cases investigated. Re- 
sults for the normalized maximum deflection ic max 
show that the value of the deflection becomes in- 
creasingly dependent upon the width-to-thickness ra- 
tio b/h as the value of the normalized end-shortening 
U increases. 

Detailed numerical results for the four theories for 
the range of b/h ratios investigated for the normal- 
ized end-shortening displacement parameter U = 14 
are presented in tables I through IV. The value of this 
applied displacement parameter equal to 14 is well 
into the postbuckling range. Results are presented in 
each table for the applied average in-plane stress re- 
sultant N xav , the maximum values of the other stress 
resultants, and the maximum transverse displace- 
ment. Table I shows that for b/h ~ 100, there is 
very little difference in the results given by the dif- 
ferent theories. Table II indicates that as the b/h 
ratio decreases, that is, as the width of the plate 
narrows, small differences between the theoretical re- 
sults occur. In tables III and IV, these differences 
become more pronounced, and although the three- 
dimensional flexibility theory does give more accu- 
rate results, the numerical procedure has difficulty 
converging for thicker plates. The omission of the 
first-order shear deformation theory results presented 
in tables III and IV and the omission of the three- 
dimensional flexibility theory results presented in ta- 


ble IV correspond to those values of b/h for which 
results were not obtained. An explanation for this 
inability to obtain results has not been determined. 

The principle of virtual work requires only that 
the geometric boundary conditions be satisfied. The 
inclusion of additional terms in the representation of 
the through-the-thickness variation of the in-plane 
displacements will lead to convergence and the sat- 
isfaction of natural boundary conditions in the limit 
if a complete set of terms is used. An alternate ap- 
proach is to use terms that satisfy the boundary con- 
ditions. A complete set of these terms also leads 
to convergence. For the present problem, the three- 
dimensional flexibility theory uses terms that do not 
satisfy the natural boundary conditions. In the 
higher order shear deformation theory, coefficients of 
u and v of the assumed displacements in equation (3) 
are chosen such that 7 xz = 0 and 7 yz — 0. The coef- 
ficients are written in terms of the existing unknowns 
u a , u a , and w° in a form which satisfies the natural 
boundary conditions at the top and bottom surfaces 
of the plate. Comparisons of results are valid whether 
or not natural boundary conditions are satisfied. 

Concluding Remarks 

This paper presents buckling and postbuckling re- 
sults for simply supported aluminum plates loaded 
in compression. The buckling results have been plot- 
ted to show the effects of varying plate width and 
thickness on the buckling stress coefficient. Buckling 
results are given for various plate length- to- width ra- 
tios. The buckling results for aluminum plates with 
width-to-thickness ratios less than 10 illustrate the 
need to include the effects of transverse shearing 
when determining the compressive buckling stress. 
Postbuckling results for plates with transverse shear- 
ing flexibility are compared with results from classi- 
cal theory for various width-to-thickness ratios. The 
plates are considered to be long with side edges sim- 
ply supported and stress free, and the plates are sub- 
jected to a longitudinal compressive end-shortening 
displacements. Characteristic curves indicating the 
average longitudinal direct stress resultant as a func- 
tion of the applied displacements are calculated based 
on four different theories: classical Von Karman the- 
ory, first-order shear deformation theory, higher or- 
der shear deformation theory, and three-dimensional 
flexibility theory. 

Present results indicate that the three- 
dimensional flexibility theory gives the lowest buck- 
ling loads for the four theories considered, and there- 
fore the most accurate results. The higher order 
shear deformation theory has fewer unknowns to de- 
termine than the three-dimensional flexibility the- 
ory but is not as accurate. The figures presented 
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show that, for postbuckling, small differences oc- 
cur in the average longitudinal direct stress resul- 
tant, in the maximum values of the other stress 
resultants, and in the maximum transverse displace- 
ments calculated when the effects of transverse shear 
flexibility are included in the various theories. Be- 
cause of small differences in the results from the var- 
ious theories, results are given in tables for the ap- 
plied edge displacements into the post buckling range. 
The principle of virtual work used in this study re- 
quires only that the geometric boundary conditions 
are satisfied. The most accurate results are obtained 
from the three-dimensional flexibility theory where 
enough terms are selected to satisfy the geometric 
boundary conditions but not the natural boundary 


conditions. In the higher order shear deformation 
theory, it is possible to select coefficients of the as- 
sumed displacements such that the transverse shear 
strains ^ xz = 0 and ^ yz = 0. These coefficients are 
written in terms of the existing unknown displace- 
ments u\ v (1 \ and w° in a form which satisfies the 
natural boundary conditions at the top and bottom 
surfaces of the plate. In this manner, results are ob- 
tained from the higher order shear deformation the- 
ory that are nearly as accurate as those obtained from 
the three-dimensional flexibility theory, but with con- 
siderably fewer unknowns. 

NASA Langley Research Center 
Hampton VA 23665-5225 
May 25, 1990 
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Appendix 

Derivation of the Governing Equations for the Higher Order Shear Deformation 
Theory and the Three-Dimensional Flexibility Theory 

Governing differential equations are derived in more detail in this appendix for the two higher order theories 
considered in this paper. 

Higher Order Shear Deformation Theory 

The displacements used in this theory are given by equations (3) as 


u{x, y, z) = u°{x, y) + <{ u a (x, y) - y 


u a (x, y) 


+ w,° ( x,y ) 


/£\2l £ \ 
\h) )h\ 


v(x,y,z) = v°(x,y) + 


f a, \ 4/i v a (x,y) ,1 / 2\ 2 1 


h 


w{x , y, z) = w°(x, y) 

Substitution of equations (3) into equations (7) gives the strain-displacement relations 
% = U°+^W,f +u “ ^ - ~ (^) 3 (u ,% +w,° xx h) 

- 3 (h) ( v ’y +w ’yy h ) 

lxy = u° +v° +w° w° + (u “ +v “ ) | ^ (^) 3 (u “ +w “ +2™ h) 


£ = W ° +-U) ° 2 f — 

£ y v >y +2 w, v +v ’V h " * • 


Ixz = ( 


-<!) 2 


lyz 


= 


ft 


1 - 


2 *y 

7T; 


The assumption of sinusoidal periodicity along the length leads to 

u° = — U — b u% (y) sin — — 

Cl A 

a a/ \ 

n = u l {y) cos — 

7TX 


w° = w°(y) sin — 


V° = vfi(y) + v$(y)cos 


27 rx 


v a = vf(y) sin 


7TX 


(Al) 


(A2) 


(A3) 
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Stresses are determined from Hooke s law according to equation (8), and stress resultant forces and moments 
are determined by the following integrals through the thickness 


N Xq 4 - N X2 cos 


27TX f 

= J- 


hj 2 
h/2 


a x dz 




27TX 

Nyo + N V2 sin — 



Gy dz 






dz 


N x , 


7 XX 

: 0 COSy 




dz 


2ixx 

N xy2 sin — 



dz 


Mxq sin ^ 


/ 


h/2 

h/2 


<T X Z dz 


(A4) 


7TX 

M T , sin — 


f h / 2 4 (z \ 3 

11 sin T = J-k/2 


nx 


r h /2 

M Un sin — = / &yZ dz 
yo A y_fc/2 


M y , sin y 




dz 




7TX 


COS 



TxyZ flZ 


Mxyi 


7TX 


COS 





dz 
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Substitution of the stresses and strains into the virtual work expression, equation (9), and performing the 
variation leads to the differential equations 


where the C{j are the stiffnesses in < 
equations, which complete the set ■ 
of the unknowns: 


'(illation (8). Using the definition ft = gives the following two differential 
equations (14 equations with 14 unknowns) without squares of derivatives 


w°f = ft 


u?/ 


|e >22 (Mj * - M Vx ) - D 22 M yi + (D22D12 - D 22 D 12 )^j 

+ (D22D12 ~ D22B l2 ) X J + < (^) 2 | } / (^22^22 - D 2 22 ) 


(A7) 


Squares of the derivatives are not retained in the computer program used. 

The boundary conditions at y = 0 and y-b that correspond to equations (10) are 

u 2 = U 1 = 0 ' 

> 

w‘{ = 0 

Myo — My 1 = ^ > 


(A8) 


Three-Dimensional Flexibility Theory 


The displacements used in this theory are given by equations (4) as 


2 „ . 7T2 

u(x, y, z) = u°(x,y) + u a {x,y)- + u' {x,y)sm — 

Z 7T2 

v(r , y, z) = v°(x, y) + v a {x , y)- + V s {x, y) sin — > 

7T2 

w(x, y, z) = w°(x , y) + w c (x , y) cos — 


(A9) 


Substitution of equations (4) into equations (7) and neglecting the nonlinear terms involving w c gives the 
strain-displacement relations 


1 n Z , . 7T2 

fx - <, +- -to “ +U “ - + ti,i Sin — 




7T2 
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e z = - R w S m T 
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7T2 
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77/2 


= t +>"°»+( r s+ <) cos Tr 


(A10) 
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The assumption of sinusoidal periodicity along the length leads to 


o tt x , ot \ • 27nr 

u = -U- + u 2 (y) sm — — 
n A 

= w f(y) cos ^ 

o o / v 7r «r 
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= v o(y) + u 2(j/) cos 
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TTX 

T 

7TZ 
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Stresses are determined from Hooke’s law according to equation (8), and stress resultant forces and 
are determined by the following integrals through the thickness 
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variation leads to the differential equations 



The stiffnesses of the plate are given by 


f h ! 2 

A-ij — I Ci j dz 
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(A14) 


where the Cjj are the stiffnesses in equation (8). The definition (3\ — w°> is not used in this theory because of 
the difference in the unknowns in the theories. Instead, the following definition of (3\ is used 


0 l 




Lyzi 
h 4 


f ^44 
\H 44 



1 + 


Ny> j f A44 
2 H 44 j \ H 44 


Hu\ 

h4 ) 


(A15) 


which results in the completely defined set of equations (16 equations with 16 unknowns) without squares of 
derivatives of the unknowns, which are not retained in the computer program used. 

The boundary conditions used at y = 0 and y — b that correspond to equation (10) are 


ul = uf = uf = 0 ' 
Nyo = N y2 = 0 
w° = w \ = 0 
My 1 — L yi = 0 


(A16) 
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Table I. Stress Resultants and Transverse Displacement for bjh — 100 at — 14 


Theory 

N x ,av 

^y,max 

N xy, max 

%iax/^ 

Classical 

7.63970 

2.69335 

1.95480 

0.15466 

First-order shear 
deformation 

7.65298 

2.70474 

1.95260 

.15465 

Higher order shear 
deformation 

7.65165 

2.73821 

1.97692 

.15571 

Three-dimensional 

flexibility 

7.61673 

2.67665 

1.96602 

.15468 


t Barred quantities are nondimensional quantities given by 


t/(.4 n -^ 2 /A 22 )6 2 

a^D^2 tt 2 


■^x,av ^ Hj _ Ny t m ax ^ -»r 

~ X!/ ' maX 


^Yry,max ^ 
V 7 ^! l£>22 tt 2 


Table II. Stress Resultants and Transverse Displacement for b/h = 20 at t {/ = 14 


Theory 

A r I ,av 

N y, max 

N xi/,max 

^max/ h 

Classical 

7.63970 

2.69335 

1.95480 

0.15466 

First-order shear 
deformation 

7.55385 

2.65594 

1.97912 

,15416 

Higher order shear 
deformation 

7.60060 

2.64025 

1.92883 

.15255 

Three-dimensional 

flexibility 

7.51680 

2.68240 

1.97278 

.15389 


t Barred quantities are nondimensional quantities given by 


F t/(ytn - /If 2 / M 22 )ft 2 
a \/^11^22 71-2 


iv. 


A r x,av b 2 

'/DuD22 71-2 


TV 


AT, 


?/, max 


j/,max 


\/DiiD 22 fl" 2 


6 2 - 

N 


xy, max 


/V h 2 

2> xy,max v 
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Table III. Stress Resultants and Transverse Displacement for b/h — 10 at U' - 14 


Theory 

-N x,av 

-V i^max 

V ;r y.max 

w ’max/^ 

Classical 

7.63970 

2.69335 

1.95480 1 

0.15466 

Higher order shear 
deformation 

7.37103 

2.53324 

1.92198 

.15039 

Three-dimensional 

flexibility 

7.20968 

2.70648 

1.99174 

.15011 


t Barred quantities are nondimensional quantities given by 


7t U(A n - A 2 2 /A 22 )b 2 

U— , — 2 ^ X,HV 

ay/ D\\D 22 7T z 


Nx,uv 

y/DuDn 


IV. 


y. inn x 


y.max b 

^11^22 7T 2 


iV 


x y.max 


jVj»y,max 

VD\\D 2 2 tt 2 


Table IV. Stress Resultants and Ttansverse 


Dis])lacement for b/h = 5 at ^ {/ — 14 


Theory 

IV*, av 

-V y.max 

^ x y.max 

U’max/ h 

Classical 

7.63970 

2.69335 

1.95480 

0.15466 

Higher order shear 
deformation 

7.26598 

i 

2.48961 

1.97121 

.14796 


t Barred quantities are nondimensional quantities given by 


tt UjAn - A\d A n)& « 

" tt 2 


N XAV b 2 
\fD\\D 2 2 7T 2 


IVy 


,rnax 


IV y.max b 2 

s/D\ 1 0 2 2 7 T 2 


IV 


xy.max 


Nj- y.max b^ 
y/DuDmir 2 
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Classical 



First-order shear 
deformation 


Higher order 
shear deformation 


3-dimensional 

flexibility 


Figure 2. Critical axial stress resultant N x ,ct versus plate length-to-width ratio a/b for different width-to- 
thickness ratios b/h and different theories for aluminum plates of finite length. 
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First-order shear deformation theory (FSDT) 
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3DFT 


HSDT 
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3DFT 
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Figure 3. Average axial stress resultant N x ^ w versus applied end-shortening U for different width-to-thickness 
ratios b/h and different theories for aluminum plates of infinite length where Ufa is finite. 
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Figure 4. Maximum transverse stress resultant N y . max versus applied end-shortening U for different width-to- 
thickness ratios b/h and different theories for aluminum plates of infinite length where U/a is finite. 
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Figure 5. Maximum shear stress resultant N xy , max versus applied end-shortening U for different width- to- 
thickness ratios b/h and different theories for aluminum plates of infinite length where U/a is finite. 
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Figure 6. Maximum out-of-plane displacement w max versus applied end-shortening U for different width-to- 
thickness ratios b/h and different theories for aluminum plates of infinite length where (7/a is finite. 
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